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ABSTRACT 

The Sculptor dwarf spheroidal galaxy has a giant branch with a significant spread 
in colour, symptomatic of an intrinsic age/metallicity spread. We present here a de- 
tailed study of the Sculptor giant branch and horizontal branch morphology, combin- 
ing new near-infrared photometry from the Cambridge InfraRed Survey Instrument 
(CIRSI), with optical data from the ESO Wide Field Imager. For a Sculptor- like 
old and generally metal-poor system, the position of Red Giant Branch (RGB) and 
Asymptotic Giant Branch (AGB) stars on the colour-magnitude diagram (CMD) is 
mainly metallicity dependent. The advantage of using optical-near infrared colours is 
that the position of the RGB locus is much more sensitive to metallicity than with 
optical colours alone. In contrast the horizontal branch (HB) morphology is strongly 
dependent on both metallicity and age. Therefore a detailed study of both the RGB 
in optical-near infread colours and the HB can help break the age-metallicity de- 
generacy. Our measured photometric width of the Sculptor giant branch corresponds 
to a range in metallicity of 0.75 dex. We detect the RGB and AGB bumps in both 
the near-infrared and the optical luminosity functions, and derive from them a mean 
metallicity of [M/H]= —1.3 ± 0.1. From isochrone fitting we derive a mean metallicity 
of [Fe/H] = —1.42 with a dispersion of 0.2 dex. These photometric estimators are for 
the first time consistent with individual metallicity measurements derived from spec- 
troscopic observations. No spatial gradient is detected in the RGB morphology within 
a radius of 13 arcminutes, twice the core radius. On the other hand, a significant 
gradient is observed in the HB morphology index, confirming the 'second parame- 
ter problem' present in this galaxy. These observations are consistent with an early 
extended period of star formation continuing in time for a few Gyr. 

Key words: galaxies: individual: Sculptor dwarf spheroidal - galaxies: stellar content 
- Local Group - infrared: stars 



1 INTRODUCTION 

The dwarf galaxies of the Local Group offer a unique op- 
portunity to quantify the evolution and interactions of low- 
mass galaxies. Dwarf spheroidals were originally thought 
to be very similar in their metallicity and star forma- 
tion histories to the galactic globular clusters, but their 
star formation history is now known to be more com- 
plex. The Sculptor dwarf spheroidal, the first dSph dis- 
covered (Shapley 1938), has a population whi ch is pre- 
dominantly old and moderately metal poor (e.g. iDa Costal 
ll984lMateoll998l.lMonkiewicz fc et alll999LlDolphml2002l) 
How ever, the presence of an extended ho rizontal branch 
(e.g. IDa Costal 11983. iMaiewski et al.l | l999ft and the detec- 
tion of associated neutral hydrogen I Carignan et alJ ll99St 



iBouchard et alJ 120031) suggest the possibility of a more 
complex star formation history. A metallicity spread has 



been inferred in Sculptor from the 


large spread of its red 


giant branch ( 


e.g. | 


Da Costa 


1984, 


ISchweitzer et alJ Il995l, 


iKaluznv et al.l 




1. IMaiewski et alJ Il999h and the pe- 



2001). Metallicity gra dients have been d i scovered in sev- 
eral dw arf s pheroidals (lHar beck et al.l200ll) . IMaiewski et all 
jl999l) and lHurlev-Keller et alJ (I1999T) showed that the red 
horizontal branch (RHB) of Sculptor is more strongly con- 
centrated towar d s the centre than is the bluer population. 
IMaiewski et alJ <ll999T) suggest this and the detection of 
two bumps in the Sculptor red giant branch are direct ev- 
idence for a bimodality in Sculptor's metallicity distribu- 



© 2004 RAS 



2 C. Babusiaux et al. 



tion. iHurlev-Keller et alJ (Il999t) find no radial gradient of 
the age or metallicity distribution w i thin S culptor, whereas 
using the same data, lHarbeck et alJ l|200lf) find a metallic- 
ity gradient. Spectroscopic studies of a few stars in the cen- 
tral regions confirmed a wide abundance range, over 1 dex 
JShetrone et al.ll20Qj) . while first results from an exten- 
sive spectroscopic study of lTblstov et alJ i2004h confirm the 
broad abundance distribution, with range -2.5<[Fe/H]<-L5, 
and show that the central regions, where the RHB stars are 
found, also contains a significant subpopulation of stars as 
metal rich as -1 dex. 

Photometric studies of the age di stributio n in Sculp- 
tor jMonkiewicz fc et alll999llDolphirj2002] and lRizzi et alJ 
120041) confirm that a significant metallicity gradient, with 
some relatively smaller age spread, are both required by the 
stellar colour-magnitude diagram, and are consistent with 
all other data. This predominately old age is puzzling, given 
that Sculptor is apparently unique among the dSph galax- 
ies in having associated HI. Sculptor is further of interest in 
that it may belong to a comm on orbital plane with several 
other dSph iLvnden-Belllll976ft . 

In this paper we present the first near-infrared pho- 
tometric study of Sculptor. The combination of our near- 
infrared data, obtained with the Cambridge InfraRed Sur- 
vey Instrument (CIRSI) on the Las Campanas 2.5m duPont 
telescope, and optical data from the ESO 2.2m Wide Field 
Imager archive, has allowed us to undertake a broader wave- 
band study of the metallicity spread in Sculptor. Observa- 
tions and data reduction procedures are presented in Section 
2. In Section 3 colour-magnitude diagrams of Sculptor are 
compared with theoretical isochrones and globular cluster 
data. In section 4 we present the detection of the RGB and 
AGB bumps. The metallicity distribution function is derived 
in section 5. Section 6 is devoted to the study of the varia- 
tion of RGB and HB morphologies with radius. We conclude 
with a discussion of the main results. 



2 THE DATA 

Photometric studies provide a valuable tool to determine in- 
ternal metallicity spreads in dSph galaxies. Although indi- 
vidual spectroscopic measurements are more precise, as yet 
the numbers of stars with direct spectra is small. Most such 
photometric studies use optical photometry. The combina- 
tion of optical and near-infrared data provides substantially 
enhanced information, reducing the effect of photometric er- 
rors, allowing colour-colour selection of sources, and provid- 
ing in particular the V-K colour, which is a good indicator 
of the stellar effective temperature. 

We have obtained wide-area near-infrared J and K-band 
photometry of the Sculptor dSph galaxy, complementing 
this with optical V and I-band data from the ESO archive. 

Figure presents the Sculptor area observed, while ta- 
bic Q summarises the observations. 

2.1 The near-infrared data 

Near-infrared observations in the J(1.25/xm) and 
-R" s (2.15/im) bands were made with the Cambridge Infrared 
Survey Instrument feeckett et al Jl997llMackav et alJl200CD 
on the du Pont 2.5m at Las Campanas Observatory. CIRSI 
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Figure 1. Field of view of the CIRSI and ESO-WFI data. The 
contours are from Irwin & Hatzidimitriou (1995). Sculptor's core 
and tidal radii are respectively 5.8 and 75 arcminutcs. 

is a mosaic imager consisting of four Rockwell lk x lk 
detectors. The pixel scale is 0.2 arcsec pixel - . The gaps 
between the detectors being comparable to the detector 
size, four dither sets are needed to fill the mosaic image, 
leading to a field of view of about 13 x 13 arcmin 2 . For each 
dither set, 9 dither frames are taken with offsets of about 
10 arcsec. Each dither frame is observed in 5 loops of 20 
seconds exposure time each. The total exposure time per 
mosaic is then 900s. Two mosaics were taken at the centre 
of the Sculptor dwarf spheroidal, as pictured in figure Q 

2.1.1 CIRSI data reduction 

The data reduction was made using an updated ver- 
sion of the InfraRed Data Reduction ( I RDR) software 
package, first developed by ISabbev et al.l l)200ll) . A sum- 
mary of the full process is given here. The updated ver- 
sion of IRDR with its documentation are available at 
|http:/ /www.ast. cam. ac.uk /~optics/cirsi/software A fuller 
description is provided bv lBabusiaux fc Gilmord (120051) . in 
their presentation of a CIRSI study of the Galactic bulge 
and bar. 

First each image is corrected for non-linearity. The dark 
current for the relevant exposure time is then subtracted. 
The data are flatfield corrected using lamp-on domeflats sub- 
tracted with lamp-off domeflats and normalized to the first 
detector sensitivity. These flatfields are also used to detect 
bad pixels and create weight maps, used during the dither 
frame coaddition. 

The sky is subtracted in two passes. A first pass sky im- 
age is made by median combining the nearest loop-combined 
frames of the dither set. After a first dither frame coaddi- 
tion, objec t masks are produced using SExtractor source 
extraction jBertin fc Arnoutsll99rf) . A masked frame is cre- 
ated from this source-detection list, with an enlarged area 
around each detected source being used. This object-masked 
frame is used to make a second pass sky subtraction on each 
loop. 

Spatial offsets between loop-combined dither frames are 
computed by cross-correlating object pixels mapped by SEx- 
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Table 1. Log of the Sculptor observations. The exposure time format (d x n x exp) lists 'n', the number of short integrations on a 
specific centre ('loops'), and 'd', the number of slightly offset ('dither') pointings. The seeing conditions, measured from the images, are 
indicated by the PSF FWHM (in arcsec) with the corresponding range during the mosaic observations with CIRSI. 



Instrument 


Field 




Filter 


Date 


Exposure (sees) 


seeing (") 


CIRSI 


Sculptor- 


-West 


J 


2001-09-04 


9 x 5 x 20 


1.02-1.4 








K s 


2001-09-03 


9 x 3 x 45 


0.8-0.9 




Sculptor- 


-East 


J 


2001-09-04 


9 x 5 x 20 


1.14-1.4 








K s 


2001-09-03 


9 x 5 x 20 


0.96-1.3 


ESQ WFI 


Sculptor 




V 


1999-07-22 


3 x 300, dit-3 


1.00 








1 


1999-07-22 


3 x 300 


0.91 



tractor. Dither frames are then coadded using a weighted 
bi-linear interpolation, excluding bad pixels. 

Finally, the astrometry is calibrated by correlating the 
SExtractor's object catalogue with the USNO-A2 catalogue. 

2.1.2 CIRSI photometry 

Once the data are reduced, we use the IRAF photometry 
routines. Source are detected using the IRAF daofind pro- 
cedure, with a significance threshold set at five-sigma. Aper- 
ture photometry is then obtained using the IRAF phot task. 
The aperture radius is assigned for each dither frame to the 
measured PSF FWHM. Observations of standard stars from 
iPersson et alJ <ll998l) were obtained each night, and were 
used to derive the magnitude zero-point of each night. Stan- 
dard star photometry used an aperture photometry radius of 
20 pixels, equivalent to a diameter of 8 arcsec. The internal 
zero point dispersion derived from multiple observations of 
the standard stars during a night is 0.013 mag in J and 0.008 
mag in K. The instrumental magnitudes derived using the 
psf-fwhm aperture are correct ed for apertu re effects using 
the curve-of-growth method iStetsorJ ll990l. implemented 
with the IRAF mkapfile task applied to selected bright 
isolated stars. Air mass corrections of fc j=0.1 and fcx=0.08 
mag per air-mass iPersson et alj H"998.f are applied to the 
photometry. Finally, images detected near the borders of 
the images are eliminated, so that only stars observed in all 
the dither frames are kept. False detections located in the 
wings of highly saturated stars are manually deleted. 

The photometric calibration was checked by correlating 
the brightest stars with the 2MASS catalogue. The CIRSI 
K s photometry is consistent with the 2MASS photomet- 
ric system. However, an offset in the J-band photometry 
°f J(cirsi) — J(2mass) = —0.042 ± 0.005 is observed. This 
is due to the 2MASS J-band filter being more extended into 
the atm ospheric water ab sorption features at around 1.1 and 
1.4 ^m llCarpenterll200ll) than is the filter used for CIRSI. 

The 5-ct magnitude completeness limits are about J~20 
and K s ~18.8. The photometric errors are about 0.025 mag 
for J<17, 0.08 mag at J=20, 0.022 mag for K a <16, and 0.08 
mag at K„ = 18.8. 

2.2 The optical data 

The optical data were obtained from the ESO 2.2m telescope 
Wide Field Imager archive. They consist of 3x300s dithered 
exposures in V and I. Each of the eight 4kx2k CCDs in the 



ESO 2.2m WFI covers around 8 arcmin x 16 arcmin of sky 
at a sampling of 0.238 arcsec per pixel, comparable with 
the CIRSI data. The total field of view of the WFI is about 
33 arcmin x 33 arcmin, including small gaps of around 10 
arcsec between CCDs. This field of view entirely covers the 
observed CIRSI fields (figure 0. 

The WFI data were processed using a variant of 
the st andard optical pipeline described by llrwin fc Lewisl 
teOOll) . After trimming, bias-correcting, natfielding and gain 
normalisation, the I-band exposures were additionally de- 
fringed using a I-band fringe frame derived from the 3 I-band 
science images. To generate the fringe frame, the 3 dithered 
I-band exposures were object masked, combined with rejec- 
tion to remove residual artifacts, and then further filtered 
to improve the local signal-to-noise ratio. 

Object catalogues were generated for each individual 
processed science frame and used to define, and refine, the 
World Coordinate System (WCS) for each frame, by com- 
parison with the online APM plate catalogues. After updat- 
ing the 2D images with the derived WCS, the 3 V and 3 I 
frames of Sculptor were stacked with cosmic ray rejection, 
using the WCS for coalignment, and confidence maps, de- 
rived from the flatfield frames and bad column lists for each 
CCD, to aid in rejection. 

Final object catalogues were then derived from the 
stacked frames, and used both to update the WCS and to 
provide morphological classification information for each de- 
tected object. Object detection proceeds via a search for con- 
tiguous pixels at a threshold above the local sky; following 
detection a series of object parameters are derived. These 
latter are used to generate position, flux and shape infor- 
matio n for use in later processing stages <lrrwinlll985l llrwin] 
Il997l) . The basic photometric measurement used is an aper- 
ture flux estimate with radius set to the average FWHM in 
the stacked frames. Additionally, all detected images have 
a series of fixed aperture measures produced (scaled to the 
basic aperture size used) in order to automatically derive 
aperture corrections for stellar images for each CCD. 

Archive observations of iLandoTtl lll992l) standard fields 
define native system zero-points in each passband using the 
colour equations for WFI available on the ESO web site. The 
gain-correction in the previous stage ensures that all CCDs 
are on the same internal system, normalised to CCD1. How- 
ever scattered light leads to a v ariation of the photometri c 
zero points across the mosaic dManfroid fc Selmanl I2OO1I . 
From the standard field observations, a correction term of 
1.5r 2 , with r being the distance in degrees from the optical 
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Figure 2. The V-J vs J-K s colour-colour diagram of all point 
sources detected in V,I,J and K s . The solid line shows the quasar 
selection criterion from the KX method (Warren ct al. 2000). 



axis, was found to be a good approximation of the effect of 
the scattered light. Overall, this provides a zero-point cali- 
bration with 1-2% accuracy. 

The 5-<t magnitude completeness limits are about 
V~24, I~22.5, so that the censorship on the data is due 
only to the C1RS1 J and K s photometric limits. The errors 
are smaller than 0.01 mag in V and 0.02 in I for magnitudes 
brighter than V~21.5 and I~19.5. 




3 PHOTOMETRY OF THE SCULPTOR DSPH 
GALAXY 

Our next task is to generate a list of point-sources which are 
stellar members of the Sculptor dSph galaxy. Extended ob- 
jects are readily eliminated from further consideration using 
the morphological flags from the optical pipeline. 



3.1 Selection of Sculptor stars 

The VJK colour-colour diagram, figure [5] allows the detec- 
tion of three other point-source populations unrelated to 
the Sculptor dSph galaxy. Red quasars are clearly seen with 
colours which are too red to correspond to any star. All 
sources with J-K s >1 were eliminated. A stream of stars 
with J-K s about 0.8 mag and V-J colour redder than 3 mag 
can be seen. They are foreground low-mass dwarfs (e.g. 
lLeggetd|l992T) : these can be eliminated, without discarding 
any Sculptor stars, by excluding all stars redder than the 
colour of the tip of the Sculptor giant branch in all rele- 
vant colour-colour CMDs (e.g. VT>1.9). A number of prob- 
able blue quasars can be seen bluer in V-J than the main 
stellar locus of Sculptor and foreground stars. These are se- 
lected and eliminate d according to the KX method criterion 
(War ren et alJl200Ci) . which is shown as the solid line in fig- 
ure |21 Remaining foreground Galactic stars are minimised 
by our selection inside each colour-magnitude diagram of 
the location of Sculptor member stars, as described further 
below. 



Figure 3. Theoretical isochrones of the Padova group (a) AGB 
and RGB isochrones for an age of 14 Gyr and for metallicitics 
[Fe/H]=-2. 3,-1. 7,-1. 3,-0.7 dex from left to right. The data points 
are our photometry, (b) The RGB and main sequence turn-off for 
ages 14 and 8 Gyr and for the same four metallicities as in the 
top panel. 



3.2 The Sculptor Colour-Magnitude diagrams 

It is apparent from figure that there is a significant real 
width to the RGB of Sculptor in the (V-K a ,K s ) colour- 
magnitude diagram, confirming several previous studies of 
Sculptor CMD. The RGB photometric width is about A(V- 
K s ) = 0.3 at magnitude K s = 16, where the photometric mea- 
surement errors are 0.023 mag. 

To determine the origin of this dispersion, theoreti- 
cal isochrones from the Padova group, given in the ESO- 
WFI and 2MASS phot ometric systems by respectively 
iGirardi et all <2002ft and iBonatto et alJ <2004l . have been 
overlaid on the Sculptor's CMDs, using a distance modu- 
lus (m -M) = 19.54 and an extinction E(B-V)=0.02 (iMated 
1998). The exti nction is derived in th e different photometric 
bands using the lCardelli et alJ 1119891) extinction curve. 

Figure shows theoretical RGB isochrones for metal- 
licities Z=0.0001, 0.0004, 0.001 and 0.004 with solar mix- 
ture, and ages 8 and 14 Gyr. It can be seen that the RGB 
is much more sensitive to metallicity than to age. The main 
se quence turn-off, which is m ore sensitive to age, was used 
bv lMonkiewicz fc et al] il999|) to der i ve an age of 15±2 Gyr 
for Sc ulptor. The results of iDolphinl (1200211 and iRizzi et all 
(2004) confirm the predominance of old stars in Sculptor. 
We therefore adopted isochrones of age 14 Gyr for figure^. 
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Figure 4. Sculptor photometric data from this study (points), 
together with the fiducial RGB lines for the sample of Galac- 
tic globular clusters from Ferraro et al. (2000), with their global 
mctallicities as defined by Ferraro et al. (1999) indicated. 



These isochrones clearly confirm the presence of a metallic- 
ity spread within the Sculptor RGB stars. 

The colour magnitude diagrams can also be com- 
pared direct ly to globular c l uster observations. IFerraro et alJ 
l|200ffl and ISaviane et alJ l)2000F) provide fiducial lines for 
the RGB of Galactic globular clusters for a wide range of 
metallicity, in the (V,J,K) and (V,I) photome t ric sy stems. 
Figure 2] shows the fiducials of IFerraro et al] i2000T) glob- 
ular clusters on the (V-K a ,K s ) colour-magnitude diagram 
of Sculptor. The transformation from absolute to relative 
magnitudes is the same as the one applied for the theoreti- 
cal isochrones. However the V and K magnitudes are not on 
the same photometric system as ours, leading to expected 
photometric differences of the order of 0.1 mag. Considering 
that globular clusters tend to have alpha-enhanced element 
ratios, whereas the Sculptor stars do not, we indicate their 
global metallicities as defined by Ferraro et al. (1999). The 
definition of this global metallicity scale and how it can be 
translated into [Fe/H] for Sculptor is discussed in section 4. 
Here again a metallicity spread is confirmed as being con- 
sistent with the width of the Sculptor RGB. Figure [I] illus- 
trates that all stars in Sculptor are more metal-poor than 
[M/H]= -1.0. 

The spread of the RGB in the (J-K S ,K S ) and (V-I,I) 
CMDs is smaller and more sensiti ve to photome t ric err ors 
than in (V-K S ,K S ). According to ISaviane et al] (l2000l) . a 
variation of metallicity from —2.0 to —1.5 dex results in a dif- 
ference in V-I of 0.04 mag at I=— 2 (one magn itude brighter 
than the RGB bump), while according to IFerraro et al] 
(2000), the same variation of metallicity at K=— 3 results 
in a variation of 0.2 mag in V-K. Considering the relative 
photometric errors in V, I and K s , this means that V-K s is 
1.7 times as sensitive to metallicity as is V-I. We will there- 
fore use preferentially the (V-K s ,K a ) CMD in the following 
to derive photometric metallicity indicators. 




4 THE RGB AND AGB BUMPS 

Local maxima are observed in the luminosity function of 
the giant branch of old metal-poor stellar populations. These 



Figure 5. Luminosity Function (LF) and cumulative LF (a) in 
the K s band for all stars, (b) in the V band separately for stars 
redder and bluer than the giant branch mean fiducial. 



RGB and AGB bumps are well known features of the colour- 
magnitude diagrams of globular clusters (e.g. IFerraro et al] 
1999). Those two bumps are detected in all our visible and 
near- infrared bands (table as illustrated in figure [5] for 
the V and K s bands. To allow us to compare the absolute 
magnitudes of these features in Sculptor with other stud- 
ies, the ESO WFI V and I photometry was converted into 
the st andard Johnson photometry using the iGirardi et al] 
(2002) isochrones in those two filter systems. At the loca- 
tion of the bumps, (V — I)j = 0.96, a given simulated star 
of the isochrones present the colours Vj — Vwfi = —0.06 
and I j — Iwfi = 0.1. As previously, the conversion to abso- 
lute magnitudes assumes a distance modulus of 19.54 mag 
and E(B-V)=0.02. Errors in the determination of the bump 
loca tion are of abou t 0.1 m ag. 

iMaiewski et al] (Il999l) also detected those two bumps 
within the central 10', but associated the second one with 
the RGB bump of a more metal-poor population. Indeed, 
the AGB bump of a population of metallicity [Fe/H]~ —1.5 
is located at the same position on the CMD as is the RGB 
bump of a population of metallicity [Fe/H]~ —2. The V 
magnitude of the second bump is consis tent with the value o f 
Mv(AGBbump) = -0.3 ±0.1 used bv IFerraro et al] <199Sft . 
Its clear detection can be explained by the fact that the lu- 
minosity level of the AGB bump stays fairly con stant with 
the cluster metallicity (e.g. lCastellani et alll99ll) . The AGB 
bump being always bluer than the RGB, it explains the de- 
te ction of this second bu mp on the blue side of the RGB 
bv IMaiewski et all 1^.999). Using, as they did, a division of 



© 2004 RAS, MNRAS OOO.EHS] 



6 C. Babusiaux et al. 

Table 2. Magnitudes of the RGB and AGB bumps. See the text for a definition of the photometric system. 



m v W Fi m hvFi m -r m i< 3 M Vj M Z j Mj M Ks 
RGB bump 19.92 18.83 18.23 17.61 0.26 -0.64 -1.33 -1.94 
AGB bump 19.40 18.38 17.75 17.13 -0.26 -1.09 -1.81 -2.42 



the giant branch into a red and blue part (figure HJ, we not 
only find as expected the second bump on the blue side of 
the giant branch, but also the RGB bump shifted by about 
0.2 magnitudes (figure^), which is consistent with a vari- 
ation in metallicity. If the second bump is to be due to a 
distinct metal-poor population, being clearly detected at all 
wavelengths, this second population should also show a clear 
imprint on the tip of the red giant branch. No such distinct 
second RGB can be detected (figures OK. and0J. We then 
conclude that the second bump is the AGB bump. 

The V magnitudes of the HB and th e RGB bump 
are g ood indicators of the metallicity (e.g. iFerraro et alJ 
1999). The peculiar shape of Sculptor's HB however makes 
the determination of its mean V magnitude unreliable (fig- 
ure [7|. The V magnitude of the RGB bump leads to a 
global me tallicity of [M/H1= -1.30±0.12 from the cali- 
bration of IFerraro et alJ <ll999l) . It s K a magni t ude l eads to 
[M/H]= -1.39 ± 0.14 according to lCho fe Leel <2002h . while 
it leads to [M/H1= — 1.19 ± 0.12 from the new calibration 
of lValenti et al. | J2004T) . Those cali brations were made using 
the relation of Isalariset alJ lll993l) : 

[M/H] = [Fe/H] + log(0.638 * 10 [a/Fe] + 0.362) (1) 

Sculptor does n ot seem to present a s trong enhancement of 
alpha-elements: IShetrone et al] i2003fl measured for 5 stars 
the followi ng values for fq/Fel: 0.18, 0.13,0.13,-0.01 and 0.23 
(table 2 of iTolstov et"al1l2003l) . For [a/Fe]=0.13, [M/H]=- 
1.3 can be translated to [Fe/H]=-1.4 by equation^ 

It should be stressed that the indicated errors do not 
take into account the uncert ainty in the d istance modulus, 
estimated to be 0.08 mag in iMateol (I1998T) . Another metal- 
licity indicator, independent of the distance modulus and 
of zero point calibration errors, is the differ ence between 
the A GB and RGB bump luminosities. From IFerraro et al] 

(1999), with WagbI^ZI = °- 52 ± °- 14 > we can derive 
[M/H]=-1.4±0.2. 

Those metallicity indicators agree with the previ- 
ous comparisons of the RGB morphology with theoretical 
isochrones and globular clusters (figure |HK and figure 0J. 



5 THE METALLICITY DISTRIBUTION 
FUNCTION 

The mean fiducial of the Sculptor red giant branch was com- 
puted through a least squares fit to a second order polyno- 
mial on the (V-K S ,K S ) CMD. Horizontal branch stars and 
foreground stars have been eliminated by selecting only the 
stars 0.2 magnitudes away from the mean fiducial (figure^. 
Only the giant branch stars brighter than K s = 18.7 and up 
to the RGB tip (K 3 >13.8) have been selected. As the AGB 
stars occupy the same location on the CMD as the most 
metal poor RGB population, no attempt to discard AGB 
stars was made. 

We derive the metallicity distribution of those selected 
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Figure 6. Least square fit of a second order polynomial to the 
Sculptor giant branch. The dashed curves frames the stars se- 
lected for the giant branch studies. 
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Figure 7. Sculptor (V-I,V) CMD. The boxes indicate how the 
HB index was computed, selecting the B,V and R stars. 

Sculptor giant branch stars using the Padova isochrones de- 
scribed previously. Each isochrone of figure EJi can be ap- 
proximated by a second order polynomial: 

V - K s = a + ai K s + a 2 K s 2 (2) 

A second order polynomial regression of those coefficients as 
a function of the metallicity of the isochrones is obtained: 

«i = G(i,o) + a(i,x) [Fe/ H] + a {lj2 ) [Fe/ H] 2 (3) 

By inverting equation each point in the (V-K S ,K S ) CMD 
can then be assigned an estimate of its metallicity. Taking 
into account only the uncertainty of the accuracy on the 
polynomical regression of equations [5] and [3] and the photo- 
metric errors, the typical uncertainty in the resulting mea- 
surement of the metallicity of a individual star is smaller 
than 0.04 dex. 
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Figure 8. Photometric metallicity distribution function as de- 
rived using theoretical isochrones. 
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Figure 9. Mean metallicity, as estimated in section 5, as a func- 
tion of radius. The vertical line show the measured dispersion 
around this value. 



The resulting metallicity distribution function is repre- 
sented in figure[H] The secondary peak at [Fe/H]= —2.2 is 
an artefact due to AGB stars and should be ignored. The 
mean metallicity obtained is [Fe/H]= —1.42 with a disper- 
sion of 0.2 dex. This mean metallicity is in agreement with 
the value obtained from the RGB bump. A comparison of 
those metallicity estimates with the literature is given in the 
discussion section. 



6 METALLICITY GRADIENT INDICATORS 

6.1 The Giant Branch 

Since the RGB in our photometric system provides a good 
indicator of metallicity, we studied its variation with radius 
as a test of a possible metallicity gradient in Sculptor. 

The metallicity indicators derived in the previous 
section have been studied as a function of radius. The 
central 10 arcminutes of S culptor has zero ellipticity 
lllrwin fc Hatzidimitrioujl 995). so as our data do not extend 
beyond a radius of 15 arcminutes, we use a simple circular 
annulus. Figure [§] show that no metallicity gradient is de- 
tected within our data. This gives an upper limit of 0.03 dex 
for the metallicity gradient within twice the core radius of 
Sculptor. 

6.2 The Horizontal Branch 

The horizontal branch of Sculptor can be clearly divided 
into a blue (B) and a red (R) part, lying on either side of 
the instability strip (V). The ratio of the number counts 
of those different parts, quantified by the HB index (B- 
R)/(B+V+R), is dependent on the metallicity. But there 
is a well known 'second p arameter problem', which could be 
age Ce.g. lLee et ai]|l994h . 

The HB index was computed from the (V-I,V) CMD, 
as illustrated in figure |7| for the same radial annuli as used 
for the RGB study. Both the full ESO WFI field of view 
data and the sub-area in common with the CIRSI field of 
view are presented in figure lTUl A K-S test for the hypothesis 
that the red and blue horizontal branch stars have the same 
radial distribution gives a significance level of 10 _9 %. This 




radius (arcminutes) 



Figure 10. The HB index against radius for all the ESO-WFI 
data and stars only within the CIRSI field of view (figure IT! . 



confir ms t he HB gradient detect ed by iHurlev-Keller et alJ 
<1999l) a ndlMaiewski et all <1999h . 

The lLee^^dT(^994Ti models give theoretical isochrones 
that show, for an HB index between -0.5 and 0.5 and a 
given age, a linear relation between the HB index and the 
metallicity: 

[Fe/H] ~ -0.34 * HBindex + cte (4) 

The observed gradient in HBindex of about 0.5 then corre- 
sponds to a gradient in metallicity of 0.17 dex. Considering 
the upper limit of 0.04 dex for a metallicity gradient de- 
rived previously from the RGB morphology, an age gradient 
is required to explain the observed HB gradient. As always 
when discussing HB morphology however, one must recall 
that the 'second parameter problem' is not yet solved and 
that another parameter may influence the HB morphology. 

The simplest conclusion is that a small gradient in mean 
age is apparent in Sculptor and that an eventual small metal- 
licity gradient associated with it would have a too small ef- 
fect on the RGB compared to the large abundance dispersion 
to be detected. 
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7 CONCLUSIONS AND DISCUSSION 

The combination of near-infrared photometry from CIRSI 
with optical data from the ESO WFI, allowed a de- 
tailed study of the Sculptor dwarf spheroidal giant branch 
morphology. We confirm that the broad giant branch of 
Sculptor demonstrates an intrinsic metallicity spread (e.g . 
Da Costalll984l ISchweitzer et al.lll99l iKaluznv etUll 19951 
M^ewskiet^]^^99Tr From our photometric study we 
quantify this spread into a metallicity range of A[Fe/H] = 
0.75 dex. The RGB and AGB bumps are detected in all 
the optical and near-infrared luminosity functions, excluding 
the substantial metal-poor cont ribution to Sculp t or's m etal- 
licity distribution proposed bv iMaiewski et ail <ll999T) . We 
derive a mean metallicity within two core radii in Sculp- 
tor of about [Fe/H]= -1.4 from both the RGB and AGB 
bumps magnitudes and isochrones fitting. Our mean metal- 
licity a nd the me t allicit y range are higher than those de- 
ri ved in |Da^jostal 1^84) from photometry of the RGB and 
m IKaluznv et all lll995l) from RR Lyrae stars: these results 
were summarized in the iMateol <ll998l) review as a mean of 
[Fe/H]= — 1.8±0.1 with a spread of 0.3 dex. However our 
photometry is in agreement with t he metal l icity e stimations 
from Sculptor RR Lyrae stars of [ Covacs|_j20pJ) and C all 
triplet observations of 37 stars of Tolstoy et al] ]200l]) . It 
is in exc ellent agreement with the very recent spectroscopic 
survey of lTolstov et a l. ( 2004), whose derived metallicity dis- 
tribution inside two core radii indicates a mean metallicity 
of —1.4 dex, and metallicity range of about 1 dex. Their 
data show the population structure is complex, with the 
more metal-poor part of the distribution function becoming 
dominant at radii beyond those we have studied here. 

We do not detect a gradient in the RGB morphology 
within a radius of 13', 2 .2 tim es th e Sculptor core radiu s. Al- 
though iH^beckeTal] ]200l]) and lTolstov et al] J2004I) find 
a metallicity gradient in Sculptor, ou r result is in agr e ement 
with their data. Indeed figure 6 of lHarbeck et all (1200 ll) 
shows that the radial distribution of blue and red RGB 
stars b egins to differ only after 13'. Figure 3 of lTolstov et al] 
i2004T) . based on spectroscopic data, indicates that a metal- 
licity gradient is indeed visible only beyond this radius. On 
the other hand, we do detect at high significance a gradi- 
ent in the h orizontal branch (HB) mo rpho logy, confirming' 
the re sults of lHurlev-Keller et al] lll999l) and lMaiewski et al] 
( 1999). As this cannot be explained by metallicity, the most 
likely second parameter could be age. iHurlev-Keller et al] 
( 1999) did not find evidence for a gradient in the main- 
sequence turn off, leading to an upper limit of a 2 Gyr 
va riation at consta nt metallicity. According to the models 
of iLee et al.l (1 19941) a small variation in age of even a few 



Gyr leads to a strong variation in the HB morphology. Age 
could then still be the second parameter. Moreover, [Fe/H] 
and age are not the only variables which can affect the HB 
morphology, but also other element a b undances, in partic- 
ular t he [O/Fe] ratio <Lee et al.lll994h . IHurlev-Keller et al] 
( 1999) detected another population presenting a gradient in 
Sculptor: a 'spur' of stars extending ~0.7 mag above the 
old main sequence turn off. They conclude that it cannot 
be explained by the presence of an intermediate age popu- 
lation, preferring the interpretation of the spur as a binary 
sequence, and speculate that it could be related to the vari- 
ation of the HB morphology. 



No significant intermediate age population has been 
found in Sculptor, excluding star formation within the last 
5 Gyr. The large metallicity spread requires extended star 
formation, while the evidence of low alpha-element enhance- 
ment implies extended star formation and self-enrichment 
over a period of at least 1 Gyr: quite long enough to affect the 
HB morphology. The HB morphology gradient implies that 
the most recent star formation episode occurred in the centre 
of the galaxy, consistent with naive expectation that gas is 
more easily retained deeper in the galaxy potential well. In- 
deed in most dwarf galaxies observed with sufficiently deep 
wide-field imaging, populations gradients have been found 
with the younger stars b eing more centrally concentrated 
(e.g. ISaviane et al.l l200ll and references therein). Our lack 
of detection of a metallicity gradient may be explained by 
the age- metallicity degeneracy that would hide a small age 
and metallicity difference and by the stronger dependence 
of the HB on other parameters such as the oxygen abun- 
dance, stellar rotation or b inar ism, as already suggest ed by 
IHurlev-Keller et al] ]l999l) and lMaiewski et al.l ]l999l) . 

All this could be consistent for Sculptor with a single 
period of sta r formation extende d in time for of the order 
of a few Gvr. lTolstov et all ]2003h conclude that their study 
of the element abundances are consistent with a closed-box 
chemical evolution scenario. The small dynamical mass of 
dwarf spheroidals such as Sculptor means that their binding 
energy is small compared to the energy released by several 
supernovae, which leads the high metallicity spread and rel- 
atively high mean metallicity derived for Sculptor puzzling: 
how did the gas stay bound long enough to have an extended 
star formation and gas enrichment? The star formation rate 
should be low to allow the chemical enrichment to proceed 
gradually. Hydrod ynamical simulations are trying to answ er 
this question (e.g. ICarraro et al-llioOll ICarigi et aljl2002ft . 
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